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FIRCLA, one-loop correction to e'^e — ^ uvH and basis of Feynman 
integrals in higher dimensions 

F. Jegerlehner and O. V. Tarasov 

^Deutsches Elektronen-Synchrotron, DESY Zeuthen, Platanenallee 6, D-15738 Zeuthen, Germany 

An approach for an effective computer evaluation of one-loop multi-leg diagrams is proposed. It's main feature 
is the combined use of several systems - DIANA, FORM and MAPLE. As an application we consider the one-loop 
correction to Higgs production in e"*" e~ —^ vVH, which is important for future e"'"e~ colliders. To improve the 
stability of numerical evaluations a non-standard basis of integrals is introduced by transforming integrals to 
higher dimensions. 



1. INTRODUCTION 

Electroweak SM calculations of one-loop cor- 
rections to processes with five and more external 
legs arc quite demanding. The number of dia- 
grams to be evaluated typically is in the several 
thousands, integrands of the individual diagrams 
are in general very lengthy and in spite of large 
cancellations which take place the final result usu- 
ally comes out to be not as compact as required 
for an efficient and reliable numerical evaluation. 
At the present time there is no computer algebra 
system which allows one to perform in an efficient 
way complete calciilations of multi-leg one-loop 
corrections. For example, FORM is very efficient 
in doing Lorentz algebra but it can not handle ra- 
tios of Gram determinants depending on several 
momenta and masses. To overcome this kind of 
problems we use a combination of several systems 
and exploit their most advantageous features. We 
are testing the effectiveness of our strategy in the 
evaluation of the one-loop correction to the pro- 
cess e"*" e~ vv H . In the last part of this note 
we also will discuss improvements possible by uti- 
lizing non-standard sets of master- integrals. 

2. FIRCLA 

To achieve a better performance we use a com- 
bination of several tools like DIANA [1] based 
on C and program libraries based on FORM 
and ?iIAPLE which also act as interfaces between 
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them. Since the main tool of evaluation is based 
on recurrence relations we called this collection 
of packages FIRCLA which stands for Yeynman 
Integral Recursive CalcuLAtor. 

FIRCLA works as follows. After the process 
is specified to DIANA it invokes QGRAF [2], to 
generate all diagrams, then constructs input ex- 
pressions suitable for use by FORM and provides 
additional information like types and masses of 
particles as well as the relations between kinemat- 
ical variables. The output from DIANA is stored 
in a file which is used as an input by our FORM 
package. Additionally, DIANA may be utilized 
to produce pictures of all or particular diagrams 
as a Postscript file. 

The FORM package performs the Lorentz alge- 
bra, takes traces of the Dirac 7-matrices, trans- 
forms products of spinors times 7 matrices to a 
chosen basis of amplitudes, and utilizing the al- 
gorithm [3], reduces tensor integrals to a com- 
bination of scalar integrals, some of them with 
shifted space-time dimension. For each diagram 
the FORM package creates a file with expres- 
sion for further processing by our MAPLE pack- 
age. The expressions are written in MAPLE for- 
mat. Each of these files also provides information 
about scalar products of external momenta and 
the masses of the particles. 

The values of all scalar invariants are calcu- 
lated in the FORM package. Relations between 
the momenta carried by the lines of the diagrams 
(pi) and the external momenta (qj) are evaluated 
by DIANA after generating a diagram and then 
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transferee! via FORM to the MAPLE package. In 
this way all useful information from DIANA can 
be transferee! to the MAPLE package. The re- 
lations between momenta are needed, for exam- 
ple, to find which integrals are independent or 
equivalent. In the FORM output files, after the 
kinematic information an expression in MAPLE 
format follows. It has the form of a sum of in- 
tegrals multiplied by polynomials of scalar prod- 
ucts, masses and spinor amplitudes. Integrals are 
just the names of MAPLE procedures like: 



npoint(JJ P(j, m?)'^-' , p, data) 
j=i 



(1) 



where P{j, nij) = (fci —pj)^ — m'j+ie is an inverse 
scalar propagator, vj its power, N the number of 
different propagators, the parameter p specifies 
the shift of the space-time dimension D = d + 2p 
and data is a set of substitutions for the momenta 
and scalar invariants. 

In the MAPLE package integrals are evaluated 
separately by using the recurrence relations al- 
gorithm for the evaluation of one-loop integrals 
described in [3,4] (see also [5]). 

If the Gram determinants of an integral are dif- 
ferent from zero then a set of three relations can 
be used: the relation for removing dots (a dot 
represents one power of momentum attached to 
the line) from lines reads 

2A„.,r/^) = ^(i+<5,,) (1^) 



d-^j.,(k-i+ + i) 



i=l 



ll^^ (2) 



Another relation reduces the shift of the space- 
time dimension and the index of the j-th line si- 
multaneously 



(a,A„) + ^(a,9feA„)k- 



fe=l 



rid) 



(3) 



Integrals with shifted space-time dimension can 
be expressed in terms of integrals in generic di- 
mension by applying the formula: 

n 

(d-^.. + l)G„_i/r^) = 



2A„ + 5^(afcA„)k- 



fe=i 



(4) 



In the above formulae the shift operators j etc. 



shift the indices 



± 1, 



A„, = 



Ym Y2n 



Yin 
Y2n 



Y 



Yij = -iPi- Pjf + ml + m], 

Pj are combinations of external momenta and mj 
is the mass of the j-th line , dj = d/dm^ and 



G. 



n-l '■ 



-2" 



PlPl 
P1P2 



P1P2 
P2P2 



PlPn-1 P2Pn-l 



PlPn-1 
P2Pn-l 

Pn-lPn-1 



When one of the determinants Gn-i and A„ is 
equal to zero it is possible to express integrals 
with n lines as a combination of integrals with 
n — 1 lines. If Gn-i = then the relation 



fc=i 



(5) 



should be applied until one of the lines of the 
integral will be removed. For the integral with 
n — 1 lines the Gram determinant may be different 
from zero and therefore the relations (2), (3) and 
(4) can be used if needed. 
If A„ = one should apply the relation 
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j=l 



k=l 



until a line of the integral will be removed. As 
in the previous case if the Gram determinants 
for integrals with n — 1 lines are different from 
zero one can use relations (2), (3) and (4). Since 
the space-time dimension of n — 1-point integrals 
will be decreased by this reduction one should 
subsequently apply the relation 



lid) 



+ 7-(rf+2) 



(6) 
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in order to increase the dimension back to d = 
4- 2e. 

Our MAPLE package reads the FORM output 
diagram by diagram, evaluates each integral from 
the input expression, adds results for the inte- 
grals and simplifies the sum. For each diagram 
the result is written in a separate file. The ex- 
pression stored for each diagram is a combination 
of master-integrals and spinor-amplitudes multi- 
plied by ratios of polynomials in scalar products 
of external momenta, masses and the dimension 
d. Summing up the results from all diagrams is 
performed by a separate program. 

3. 0{a) CORRECTION TO e+e" ^ Vi^H 

In order to check the effectiveness of our 
method of evaluating Feynman diagrams we cal- 
culated the one-loop correction to Higgs produc- 
tion via 

e"^('3'i)e~('72) ^ V{q:>.)v{q^)H{qi) . 

This process will be important at future e"'"e~ 
colliders. In the tree approximation it was con- 
sidered in [6] . At one- loop order recently in [7] . 

At e+e~ linear colliders operating in the 300- 
800 GeV energy range, the main production 
mechanisms for SM-like Higgs particles are 

e+e- {Z)^Z + H 

e+er V V {WW) v v + H 

e+e" e+e"(ZZ) e+er + H 

i.e., "Higgs-strahlung" , VFT4^- fusion (see figure), 
ZZ-iusion, and "radiation off-top" , respectively. 
The first two processes 

e+ H e+ 



e 




H 



are the dominant ones, in particular at energies 
above 500 GeV. 

The Higgs-strahlung cross section for large s goes 
like ~ 1/s and dominates at low energies. In con- 
trast, the lyVF- fusion mechanism exhibits a cross 



section growing like ~ \og{s/M^) and therefore 
dominates at high energies. At -y/s ~ 500 GeV, 
the two processes have approximately the same 
cross sections. The relevance of one-loop correc- 
tions was considered in [8]. 

In order to generate all required diagrams a 
'technical SM model' with two leptons doublets 
(in order to distinguish external leptons from lep- 
tons in loops) and one quarks doublet was consid- 
ered. The contributions of the missing fermions 
have been obtained by adding up the correspond- 
ing results with the appropriate masses. 

We put TOe = but keep all other masses 
m^, mfc mw , vnz-, mt, mu different from 
zero. The results for the pentagon diagrams de- 
pend on 5 scalar invariants. All calculations were 
done for arbitrary values of the gauge parame- 
ters iz, $7- In the specified model 326 dia- 
grams contribute to the one-loop correction. Out 
of these 15 are pentagon type diagrams. 

The results have been stored for each diagram 
separately in a file. Its size is ^ 100 MB ( 
~ 5Mi? in MAPLE format). The expression for 
each diagram has the form 

where li are master integrals and Oj are 24 spinor 
amplitudes. 

After summing all diagrams we find that the 
result is gauge invariant as it must be. Each pen- 
tagon diagram separately gives a gauge invariant 
scalar 5-point integral - 15 in total. We observed 
a huge reduction in the number of scalar n-point 
(n=4,3,2) integrals Dq, Cq^Bq. For example, from 
240 integrals Co only 79 remain in the final an- 
swer. It should be mentioned that at the initial 
stage before taking into account the symmetries 
of the integrals the total amount of master in- 
tegrals Co with different values of momenta and 
masses was more than 500. A special routine has 
been added to the MAPLE package which recog- 
nizes the symmetries of integrals. By permuting 
masses and momenta into lexicographical order, 
taking into account momentum conservation, one 
obtains a set of independent standard integrals 
and hence a much more compact representation. 

The final gauge invariant result is stored in a 
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file of about 1.2 MB in MAPLE format. Thus, it 
is still rather large for further numerical evalua- 
tions. Working out a more compact representa- 
tion of the result, mainly by optimizing ratios of 
huge polynomials depending on several variables 
is one of the problems. 

To obtain the complete radiative correction to 
the observable differential cross section we still 
have to add the bremstrahlung correction to our 
result. This work is in progress. 

One observation we made is that in the final 
result the pentagon diagrams yield terms of the 
form 
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(d-4) 



r(d) 

'mil 



(7) 



In fact, second rank n-point tensor integrals pro- 
duce terms of the form (jn^qrjy-fi?..! with a coef- 
ficient which is singular at d = 4 for n = 5. In our 
case, these terms originate from 2nd rank tensor 
integrals like 



/ 



.ki.-.V) 



C1C2C3C4C5 



(8) 



The tensor integral here can be written in terms 
of scalar integrals with shifted dimension: 



/ 



d^ki- 



C1C2C3C4C5 



2i'A"' -'11111 



+ 2^11311 P.>P3i^ + 2/;ii3/P4mP4i. 



+42m {pi>P2}m/v + -^2121*1 {Pi.-Pa}^.^ 

+4im {P1)P4}mi' + Il221l{P-2,Pi}^iv 

+-fi2m {^'2,^4}^!/ + 111221 {Ps,Pi}^iv 



r(d+4) J 



where {pl,p2}^,u = Pi^iP2u + PiuP2ti- 

Different tensor structures give contributions 
to different of the abovementioned spinor ampli- 
tudes. By using the recurrence relations one gets 



it looks like we have to evaluate the e term in the 
expansion of the pentagon integrals. There are no 
such problems for the 3- and 4-point functions. 

In fact multiplying (8) by the Born term, sum- 
ming over polarizations and taking the traces re- 
moves the problematic — 4) terms. However, 
if we would attempt to calculate the amplitude 
numerically before squaring it, we would have to 
work out first it's e-expansion. A possibility to 
avoid this problem is to use a different basis of 
master-integrals in higher dimensions. 

4. MASTER-INTEGRALS IN HIGHER 
DIMENSIONS 

The idea to express one-loop tensor integrals 
in terms of integrals with shifted dimension was 
proposed in [9] and later in [5]. Recurrence rela- 
tions which allow us to reduce any one-loop inte- 
gral with shifted dimension to a combination of 
integrals in generic dimension were given in [3] 
and [4]. These were extended to massless inte- 
grals with both Gram determinants zero in [10]. 

In [11] it was discovered that one- loop inte- 
grals in higher dimensions provide better numer- 
ical stability. 

A simple formula expressing any n-point inte- 
gral in terms of integrals in higher dimensions was 
given in [4]. It reads 

jid) ^ {d-n + l)Gn-i (rf^2) y-^ i^k\^)_ (d) 
2A„ " 9A„ 



k=l 



2A„ 



The improved stability of numerical integrations 
of integrals in higher dimensions can be seen from 
their integral representation 



where 



Ad+2) 

'11111 



2A,, 



{d ~ 4)Gz 



-'11111 + / , 



{d - 4)G4 



(d) 



The integral at the left hand side is UV and IR 
finite, however, in front of the integrals at the 
right hand side the spurious pole 1/s appears. So 



J2 
h2 
h3 



1 



Jn 



-2-i„_3- • -J-i 



Pl2^1 



„2 



12 

-a;ia;2(l - xi)p{^ - xlx2{l - 2:2)^^2 



■ml) 
•clx2{l - X2)p'^ 
-xi (1 - a;i)(l - X2)pl3 + xiX2m\ 
+xi{l — X2)m% -I- (1 - x-i)m\ . 
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If we transform integrals to the dimension D = 
d+ 2n — 2 (assuming that d = 4 — 2e) then the 
expansion of 7^''+^""^) for smaU e = 2 — d/2 reads 



j{d+2n-2) = r (1 + e) 



e 



K + 0{e) 



where 



9n= [ ■ ■ [ {dx} J„hn = — 

Jo Jo (" 



+ 1)! 



1 

--T 



RI, 



{dx} Jnhn\nh„. 

If we transform the integrals to dimension D = 
d+2n- 
e reads 



d+2n—4: then the expansion of I^'^^"' at small 



j{d+2n-4) = r (1 + £) 

where 



1 



{n-iy.s 



Rl + 0{e) 



Rl= I ... I {dx) J„ln/i„. 
Jo Jo 

We see that in the integrals there are no 
polynomials in the denominator and this is why 
these integrals in higher dimensions are more suit- 
able for a direct numerical evaluation. 

As an example we give the relations for some 
integrals in the logarithmic basis. In the latter 
for n > 2 we express the integrals I^^ in terms of 
integrals : 



-52Ai234''\l3)-93Ai234'''(12), 

4Ai234Ai23Ai24Al34A234-f4''' = {d — — 3) 

^ 5123 Al23 Al24 Al34 A234.f4''^^^ 

j_A3) , .(3) ,(3) .(3) ,(2) 

"r^l234 "T J2134 "T J3124 J4123 ./l234 
,j-(2) , ^(2) .(2) .(2) „(2) 
"r/3124 "T /4123 "T /3214 "T /4213 "T /4312 



(d+2) 



9iAi234''^(23) 



where 



+ Ajj j2i3i4^n Aiii3j4 Aj2i3j4) 

^ Ajjj2i3i4AiiZ2i3 A7^72i4-^2 (^3^4) 



filhiaii ~ (^ii^iii^iaiiiiA '^)9i2i3ii^iii2i3i 

+ (d — 3) 9ji j2i3i4 Ai2i3i4 ) 

\ \ \ 7-(<i+2)/- • ■ \ 

^ ^nt2«3'^ni2j4'^ili3J4-'3 U2*3*4j- 

where the following notation is used: 



A,: 



= A, 



9ii...i„ 



Gn- 



1 ■ 



In these expression we see multiple occurrence of 
Gram determinants and their derivatives. This 
fact may lead for some problems for kinemati- 
cal regions where Gram determinants are close to 
zero. However outside these regions we expect 
good numerical stability of these integrals. In 
fact there is a correspondence between the inte- 
grals in higher dimensions and the method advo- 
cated in [12] for a direct numerical calculation of 
loop integrals. For the one-loop case, (3) and (4) 
are the explicite solutions of the relations which 
are derived in [12] by exploiting the Bernstein- 
Tkachov theorem. 
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